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Femtosecond time-resolved photoemission of photoexcited 
electrons in highly oriented pyrolytic graphite (HOPG) pro- 
vides strong evidence for anisotropies of quasiparticle (QP) 
lifetimes. Indicative of such anisotropies is a pronounced 
anomaly in the energy dependence of QP lifetimes between 
1.1 eV and 1.5 eV - the vicinity of a saddle point in the 
graphite bandstructure. This is supported by recent ab initio 
calculations and a comparison with experiments on defect- 
enriched HOPG which reveal that disorder, e.g. defects or 
phonons, increases electron energy relaxation rates. 

Studies of the ultrafast dynamics of electronic excita- 
tions in solids have challenged experiment and theory for 
years. Advances in ultrafast laser technology and the 
latest theoretical developments allow us to investigate 
electron dynamics in growing detail and further our un- 
derstanding of fundamental scattering processes in solids 
at the femtosecond time scale 

For an interacting 3D-electron gas the standard the- 
ory of e-e scattering - Landau's theory of Fermi-Liquids 
- predicts a quadratic dependence of scattering rates on 
the quasiparticle (QP) energy (E — Ep) ||- In a peri- 
odic potential, however, the electronic states are modi- 
fied with respect to those of a free electron gas and form 
Bloch states which may result in different QP lifetimes 
at the same energy if distinct fc-states within the Bril- 
louin zone are compared. Such anisotropies were indeed 
predicted by recent ab initio self-energy calculations for 
Al and Cu, which find strong variations of QP lifetimes 
for electrons in different bands ^|,|5| . Experimental verifi- 
cation of these predictions, however, remains truly chal- 
lenging. Time-resolved photoemission from different cop- 
per surfaces has revealed some dependence of the mea- 
sured electron dynamics on the crystallographic surface 
orientation, but the observed effects could not be clearly 
attributed to anisotropies of QP lifetimes |6| . Similar ex- 
periments on aluminum [Q have likewise not been able 
to identify band structure effects predicted theoretically 

Graphite, a semi- metal with layered structure, is 
expected to be an ideal candidate if band struc- 
ture effects are to be observed experimentally. The 
strongly anisotropic band-structure of graphite is ex- 
pected to furnish electron scattering processes with simi- 
lar anisotropies, leading to anomalous QP lifetimes. The 
special topology of the graphite bandstructure has indeed 



been used to explain apparent deviations of the electron 
dynamics observed on a ceasiated HOPG surface || from 
the standard predictions for a 3D-electron gas II . How- 
ever, the experiments provided no evidence for anomalies 
that might be associated with anisotropic QP lifetimes. 
Furthermore, the overlap of the energy range probed ex- 
perimentally in Ref. || and the range of validity of the 
approximations used in Ref. (9) is small and calls for 
more extended and more detailed theoretical as well as 
experimental work. 

Here we present a time-resolved photoemission study 
of the electron dynamics in HOPG, with particular focus 
on its energy dependence and the influence of disorder 
on the relaxation rates. The most striking observation 
- a pronounced anomaly in the energy dependence of 
the electron dynamics between 1.1 eV and 1.5 eV - can 
be associated with electrons near a saddle point in the 
graphite band structure at the M-point of the Brillouin 
zone. This, and a comparison with experiments on defect 
enriched HOPG provide strong evidence for anisotropies 
of QP lifetimes in graphite which have also been pre- 
dicted by a recent ab initio calculation JTT[ | . 

The HOPG sample Jl| was attached to a tantalum 
disk which could be temperature controlled from 25 K to 
1200 K. It was mounted in an UHV chamber with a base 
pressure of 2 x 10 -10 mbar. Experiments were performed 
at room temperature, unless mentioned otherwise. The 
HOPG sample was cleaved directly before being trans- 
ferred into UHV and heated repetitively to 900 K. The 
defect enriched HOPG surface was produced by Ar-ion 
sputtering (0.5 keV, 7 fiA, 300 s). 

Photoelectron spectra are obtained by means of the 
time-of-fiight technique with an energy resolution of 
lOmeV. For time resolved pump-probe measurements a 
visible pump and frequency doubled probe pulse of typi- 
cally 85 fs duration are focused on a spot of 50 fim diame- 
ter on the sample. The photon energy of the probe pulse 
was chosen to exceed the work function of the sample 
(e$ = 4.50 eV ±0.05 eV) by typically 0.15eV in order to 
probe not only low energy electron but also hole dynam- 
ics. A more detailed description of the experiment can 
be found elsewhere [ fL2| . 

We start our discussion with a brief overview of pro- 
cesses contributing to the dynamics of photoexcited car- 
riers. After optical excitation the energy absorbed by the 
electrons is redistributed via electron-electron scattering 
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(e-e), electron-phonon scattering (e-ph) and transport 
from the surface into the bulk. In graphite the trans- 
port of electrons out of the photoemission detection vol- 
ume into the bulk is hampered due to the comparatively 
weak interlayer coupling and is further suppressed by fre- 
quent stacking faults in HOPG Therefore, the dy- 
namics observed in the experiments are expected to be 
dominated by e-e and e-ph scattering processes with e- 
ph scattering becoming increasingly important closer to 
the Fermi level. 
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FIG. 1. Cross correlation trace: change of the photoemis- 
sion intensity within a certain energy interval as a function 
of the pump-probe time-delay. The pump photon energy was 
2.3 eV. The solid line is the fit to a bi-exponential decay. Fast 
and slow components are indicated by the dotted and dashed 
lines respectively. The inset shows a 2PPE-spectrum for si- 
multaneous pump-probe excitation. 

A typical photoelectron spectrum for simultaneous 
pump and probe excitation is shown in the inset of 
Fig. |l]. The electron dynamics are obtained by plotting 
the change AI of the photoemission intensity, resulting 
from excitation by the visible pump pulse, as a function of 
the pump-probe time-delay (see Fig.Q). These cross cor- 
relation traces (XCs) can clearly be characterized using 
a bi-exponential decay with a fast component decaying 
on the sub-picosecond time-scale and a slower component 
decaying on the ps time-scale. 

A detailed analysis of photoelectron spectra obtained 
at different pump-probe time-delays shows that the elec- 
tronic system approaches an internal equilibrium with 
a characteristic time constant of (250 ± 50) fs. This is 
estimated using the deviation of photoelectron spectra 
from the best fit to a Fermi-Dirac distribution Jl5| , p7t . 
This process - referred to as internal thermalization - is 
associated with the initial fast decay seen in XC-traces 
and results in an increase of the electron gas tempera- 
ture by 400 K. This is in reasonably good agreement with 



the calculated temperature rise using the electronic heat 
capacity of graphite and the estimated laser power de- 
posited within the optical penetration depth. The slow 
component of the XCs can be associated with thermal- 
ization of the electronic system with the lattice due to 
e-ph interaction ]l7|-|l9||. 

In the following we focus on the initial fast decay, asso- 
ciated with the internal thermalization of the electronic 
system. The corresponding decay rates - as obtained 
from the fast component of the bi-exponential fit - are 
plotted in Fig. || as a function of the intermediate state 
energy (E — Ep) = Eki n + e<f> — hv pro i, e . The data shown 
in Fig.|] extend to energies ~ 0.15 eV below Ep, and 
present, to the best of our knowledge, the first 2PPE 
measurements of electron and hole dynamics close to Ep. 
Note, that lifetimes of holes in the copper rf-bands have 
been deduced indirectly from hot electron lifetimes by 
Petek and coworkers JTgJ . 
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FIG. 2. Energy dependence of decay rates on the pristine 
HOPG surface. The anomaly near 1.5 eV is indicative of 
anisotropic QP lifetimes. This is supported by a recent ab 
inito calculation of QP relaxation rates for graphite which 
shows a similar anomaly in the same energy region [11]. 

Obviously the energy dependence of the decay rates 
cannot be described by a (E — £?F) n -dependence with 
ri = 2 as expected for a Fermi Liquid or n — 1 as previ- 
ously predicted for a simplified bandstructure model of 
graphite ||. The most striking deviation from a simple 
power law is the anomaly between 1.1 eV and 1.5 eV. In 
this energy range we observe a plateau-like region with 
no significant increase of the relaxation rates. The error 
bars in Fig. 2 reflect statistical errors obtained from the 
bi-exponential fit after averaging over 13 experimental 
runs. Systematic errors due to day-to-day variations are 
somewhat larger. 
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FIG. 3. (a) Graphite bandstructure along high symmetry 
lines and density of states, (b) e-h pair spectrum for a free 
electron gas and a graphene sheet (left panels) as well as 
energy-loss vs. momentum-change for electrons decaying from 
an intermediate state with E t = 1.5 eV (right panels). The 
overlap of e-h pair excitation- with energy loss-spectra is given 
by the hatched area. 

The most obvious feature in the graphite bandstruc- 
ture with which this plateau may be associated is a sad- 
dle point of the 7r*-bands at the Af-point of the graphite 
Brillouin zone (see Fig.||a)) and the resulting van Hove 
singularity at 1.5 eV in the DOS. Electrons excited to 
this region of the band-structure can only decay by a 
transfer of comparatively large momentum to secondary 
excitations {e-h pairs), in contrast to electrons decay- 
ing from other iso-energetic regions of the bandstructure. 
This is illustrated schematically in the lower right panel 
of Fig.|]b) for a simplified band structure model of a 
graphene sheet JlO[ | . The e-h pair excitation spectrum 
of graphene on the other hand (lower left panel) has no 
overlap with the energy-loss and associated momentum- 



change for electrons decaying from the vicinity of the Ap- 
point. In contrast, a much smaller mismatch exists for 
iso-energetic electrons along the T — K line. We speculate 
that this gives rise to strongly anisotropic QP lifetimes in 
graphite and ultimately leads to the observed anomaly in 
the decay rates. Note, that no comparable anisotropics 
exist for a free electron gas (upper panels in Fig^b)) or 
for lower energy excitations in graphene. 

Recent ah initio calculations of QP lifetimes in 
graphite by Spataru et al. pd[ | show a similar anomaly 
around 1.5 eV. The authors find pronounced anisotropies 
in the QP lifetimes with a distinct anomaly at 1.5 eV aris- 
ing from electrons in the 7r*-band around the M-point of 
the Brillouin zone. If averaged over the entire Brillouin 
zone - mimicking the averaging process in 2PPE due to 
scattering in the intermediate and final states - the cal- 
culated decay rates are reproduced by the dashed line in 
Fig. 2. The qualitative trends and features of both, ex- 
perimental and theoretical decay rates are evidently the 
same. 

The calculations by Spataru et al. strongly support, 
that the experimentally observed plateau around 1.5 eV 
is a direct consequence of anisotropic QP lifetimes in 
graphite. A quantitative comparison of theoretical and 
experimental results, however, appears to be difficult. 
This is most likely because theory does not account for 
finite temperatures, unavoidable lattice defects, and pos- 
sible corrections from diffusive or ballistic transport of 
excited carriers out of the detection volume; effects which 
are expected to lead to higher decay rates as observed ex- 
perimentally. 

The influence of disorder, e.g. defects or phonons, on 
the electron dynamics is expected to depend strongly on 
anisotropies of QP lifetimes. This is best understood 
if one considers that lattice distortions break the trans- 
lational symmetry of the crystal lattice and, thereby, 
lead to a coupling of different electronic states. If the 
states coupled in this manner have different lifetimes this 
should ultimately change the global electron dynamics. 
Matthicscns rule - which states that the rates of compet- 
ing decay channels simply add up - would imply that the 
global dynamics is increasingly dominated by short-lived 
states once the coupling among states in the Brillouin 
zone is turned on. 

To test this hypothesis we introduced lattice defects by 
Ar-ion sputtering. In Fig. ^ we plot the relaxation times 
of the electron population - the inverse of the rate - as 
a function of intermediate state energy. Qualitatively we 
find that the decay of the electronic population is accel- 
erated significantly irrespective of the electron energy. 

In agreement with the arguments presented above, 
these results provide further evidence for anisotropies in 
the QP lifetimes in graphite. Furthermore, they extend 
a recent study on the influence of defects on the electron 



dynamics in surface states 1 20 to the dynamics in bulk 
states. 
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FIG. 4. Influence of defects and lattice temperature on the 
electron dynamics. 

Alternatively, coupling between electronic states can 
be introduced by heating of the sample, i.e. by an in- 
crease of dynamic lattice distortions. This is most clearly 
illustrated by plotting the decay-time found at the peak 
of relaxation time curves like the one in Fig. as a func- 
tion of lattice temperature (see inset of Fig.||). Very 
much like static defects, lattice vibrations can lead to a 
coupling of electronic states and reduce the measured re- 
laxation times. We note, however, that the trend seen in 
Fig. [| may be a superposition of different effects with ad- 
ditional contributions from the temperature dependence 
of e-e scattering rates and enhanced inelastic e-ph scat- 
tering. 

These results may also shed new light on the appar- 
ent discrepancies between electron dynamics measured 
by different groups in graphite. Xu et al., for example, 
have observed somewhat shorter relaxation times on a 
cesiated graphite surface ^ whose absolute values lie in- 
between our results for the pristine and the sputtered 
surface. Furthermore, as in our experiments on the de- 
fect enriched surface, neither Xu et al. nor Ertel et al. 
pl[ find any evidence for a plateau around 1.5 eV This 
may indicate that surface preparation and defect density 
plays a more important role for measurements of bulk 
electron dynamics than previously assumed. 

In summary, we have investigated the electron dynam- 
ics in HOPG resulting from ultrafast optical excitation. 
The observed anomaly in the energy dependence of the 
electron dynamics can be linked to strongly anisotropic 
lifetimes of electrons excited to about 1.5 eV above the 
Fermi level. Electrons excited to a saddle point in the 
graphite bandstructure at the same energy are expected 
to be very long lived due to a mismatch of their en- 
ergy loss spectrum and the e-h pair excitations available 



for scattering processes. Disorder is found to accelerate 
global electron dynamics, most likely due to the coupling 
of long- to short-lived states. This provides further evi- 
dence for the anisotropy of QP lifetimes, which have also 
been predicted by a recent ab initio calculation [ p"l| . 
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